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Abstract 


Interactions between galactic cosmic, rays and matter are a prim ary focus of the ft ASA 
radiation problem. The electron rai/nct ic forces involved are for the most part, well doc- 
umented. Building on previous research, this study investigated the relative importance 
of the weak forces that occur when a cosmic ray impinges on different, types of materi- 
als. For the familiar elect mniagne.tu : rn.se. it is known that energy lost m the form of 
radiation is more sigriificant than that lost via contact collisions, the rate at which th< 
energy is lost, is also well understood. Similar results were derived for the weak force 
case. It was found that radiation is also the dominant mode of energy loss m weak 
force interactions, and that weak force effects are indeed relatively weak compared with 
electromagnetic eff rets. 

1. Introduction 

When a fast-moving charged particle passes through matter, it loses energy and slows down. 
There are two ways it can lose energy. It can either lose kinetic energy during collisions with 
other particles or lose energy in the form of radiation. During any given collision, the incident, 
particle will transfer some' of its kinetic energy to other particles and thereby undergo an accel- 
eration. It is known that whenever a charged particle is accelerated, it radiates electromagnetic 
(EM) energy. This type of radiation is referred to as bremsstrahlung (braking radiation) because 
it was first Observed in experiments where high-energy electrons were stopped in a thick metallic 
target. The theory of energy loss by moving charges is well understood for EM interactions. One 
important result from these studies is that for nonrelativistic particles, the energy lost in the form 
of bremsstrahlung is negligible compared with the kinetic energy loss, whereas for relativistic par- 
ticles, bremsstrahlung can be the dominant form of energy loss. This paper addresses the physics 
of energy loss by high-energy cosmic rays as they travel through bulk matter, such as the wall of 
a spacecraft or the Earth’s atmosphere. Thus, the bremsstrahlung processes were of central concern. 

For most macroscopic particle phenomena, such as the passage of a high-energy particle through 
bulk matter, energy transfer via. EM processes is usually the most, significant channel, but contribu- 
tions from weak force processes can sometimes play the dominant role. The purpose of this project, 
is to generalize the analysis of energy loss by charged particles via EM interactions to that via weak 
interactions. In particular, one quantity was sought after for both bremsstrahlung and collisions! 
energy loss mechanisms stopping power (SP), dE rM m/dX. The SP is the energy lost by the 
incident particle per unit path length traversed. The SP generally depends on the energy, mass and 
charge of the incident particle, and the density and charge of the target, material. 



2. Method 


2.1. Previous research 

The method used here to find the SP via weak interactions for the two types of energy-loss 
mechanisms was to follow the well-documented formalism for similar EM processes. The reason for 
this approach was that the equations for weak interactions differ from those for EM interactions by 
only two parameters: the charge Q of the particle undergoing the acceleration and the mass M of 
the particle mediating the force. Formulas wore appropriately generalized to cases where Q and M 
differed from the EM case. 


For EM interactions, Q is the familiar electric charge, which is always an integral multiple of 
the charge — c of the electron. For weak interactions, Q is known as the weak charge and is only 
slightly more complicated than the electric charge. In all cases, Q is deduced from a knowledge 
of the 4-current J 1 ' — ( p, .7 ) that describes the charge p and current ,7 densities of the particle in 
question. The correct form for J'' for any given interaction is provided by the Standard Model. For 
all elect roweak interactions of interest, the 4-current, of a point charge in motion is found to be of the 
form ,P‘(r,t) = <£[r(/)]c/'', where $[r(t)] is the usual Dirac delta function, which vanishes everywhere 
except, at f{ t ) = 0, and <f is a new quantity called the 4-cliarge. Generally <f is found to be of 
the form </'' — (p u 1 ' + where u'' and s ft are the 4- velocity and normalized 4-spin, respectively, 
of the particle and (j\ ■ and q\ an' coefficients (the vector and axial-vector charges, respectively) 
that, depend on various well-known parameters (viz, electric charge, weak isospin and weak mixing 
angle) that appear in the Standard Model. In a frame in which the particle’s velocity is v = vz and 
spin direction is s, w" = 7(1, v) and *" = 7(Au,.s), where 7 = l/y/1 - a 2 is the Lorentz factor and 
A = .s • 2 is the normalized helieitv of the particle. Note that if the direction of “spin-up” is taken to 
be + 2 , then A will be either +1 or -1, and, further, in the ultrarelativistic limit of interest (where 
v — > 1), .s'' = A»'', and t.lms ./''(r, t) - fi[r(0]('7 + A^)//''. Using this expression for the 4-current, 
a formula for the 4-potential A* = ( < k,,4), which describes the scalar 4> and v(>ctor A potentials of 
the point charge, was deduced (as a solution to the Proca wpiation). Without specifying the boson 
mass M at this point, ,4'' is related to by the following equation: 
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The priiiK'd <piant,ities in this equation label the (differential) source charge elements. By using 
the expression for ,P‘ involving the Dirac 1 delta function, the integral is easily seen to reduce to the 
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Electric E and magnc'tic B fields are found from these equations in the same way that they are in 
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the massless ease: 
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For a point charge' moving at v ~ 1 in the i din'ction, tlie only nonvanishing components of the 
potentials and fields, evaluated at impact parameter />, aie 
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where r = yjl? + (7c/) 2 . The quantity of central concern to energy transport by the fields is the 
Povnting vector, S. For massive fields traveling through a vacuum, S can be shown to be generally 
of the form S = E x B + M 2 A>A. Clearly it will always depend quadratic-ally on the charges. 
Furthermore, as seen above, the helicity A of the point charge' invariably enters into this factor. As 
the usual application of any such theory is to unpolarized beams of particles, an average overall 
possible helicity state of the charge is of interest. It turns out that the charge-related quantity 
that appears in the final set of equations describing the distribution of bosons outside the moving 

charge is the quantity q(- + <&. Q is thus defined as Q = y/q(, + <ll Q is Ilot th(> nsual < ‘ anolli< al 
charge (which happens to be simply qy), as defined via the Noether prescription, so it can only be 
referred to as some sort of (Lorentz invariant) effective charge. For electromagnetic interactions, 
qy = and q A - 0, where e = y/4na is the electromagnetic coupling constant (a 1/137 is the 

fine structure constant) and Q VM is the familiar dimensionless charge' quantum number for electro- 
magnetic interactions (e.g., (? EM = ±1 for positrons and electrons, respectively). For neutral weak 
interactions (mediated by Z bosons), q v = g z {T - 2Q EM sin 2 0 w )/2 and q A = ~9xT / 2, where 
9\y ~ 28.74° is the weak-mixing angle', T* is the third component of the vectoi of weak isospiu 
quantum numbers of the point charge (e.g., T :i = ±1/2 for neutrinos and electrons, respectively), 
and g z - c/ sin^r cos^r is the neutral weak coupling constant. And finally, for charged weak 
interactions (mediated by \\ r± bosons), qy = gw/ 2\/2 and q A = =F.9 iv/ 2\/2, respectively when* 
g w = v j sin die is the charged weak coupling constant. 

One guiding principle for determining M was conservation of energy and momentum. So M 2 = 
(AEj 2 — (A p)' 2 , where A E and A/7 are the energy and momentum transferred (mediated by the 
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boson) in the collision, serves as a defining equation for this parameter. The boson masses are 
not uniquely defined in general by this equation, however, and an additional piece' of information 
is needed. The Lorentz condition in momentum-space, p tl A> 1 = 0, yields a relation between AE 
and A/7 (here //' - (AE, A/7) is the boson’s 4-momentum) that can be used to eliminate this 
ambiguity. It was actually one of the more' interesting results of this study that AE and A/7- (the 
/-component of A/7) are found to be related to one another by wav of the above-mentioned charges: 
Alibi ~ AE(l — sj where z = q A A/(q\’V + \q A ) is a new charge parameter. After taking the 
ultrarelativistic limit, averaging overall possible helicity states of the colliding particle, and making 
other slight reparameterizations that ensure M is Lorentz invariant and the concept, of causality is 
preserved (i.e., the boson travels subluminally), an equation for M was arrived at that is unique to 
each of the tlnce types of interactions of interest.. The final mass scheme is as follows: 


M = M f 





(5) 


A new charge parameter, appears here. It is very much like e defined above, and. in fact , is 
derived from s. For the photon, // is found to be 7/ = q A /(q A ~q{) = 0, so that M = (), as one might 
expect.. For the Z boson, // = q A /(q A - qv) = T :i (T :i - Q EM sin 2 0 H -)/2, and 7/ = 2 q A /(q A T q v ) = 1 
for the U'= boson, respectively. Interestingly, the masses of the W and Z bosons do not equal 
the' masses predicted by the Standard Model (80.419 and 91.1882 GeV, respectively), and in fact 
can be orders of magnitude different, yet the distribution functions for the bosons using this mass 
scheme agree exactly with the well-known distribution functions found in the parton theory for 
electroweak interactions. It is also noteworthy to point out that the ratio of the masses of the Z 
and YV bosons here is roughly the same as the ratio found in the Standard Model. For example, if 
the colliding particles are either protons or electrons, this ratio is 0.9300 lien', and cosd w = 0.87G8 
in the Standard Model. The other parameters found in the above equation are as follows. A It is 
the mass of the particle (/ for fermion) that emits the boson and x = A E/E f (E f is the energy 
of the particle) is the Feynman scaling variable, which is bounded (by conservation of energy and 
causality) by the inequality string qv/\(yv) 2 + r/] < ;r < 1 - I/7. I11 summary, a boson mass 
scheme was derived from the basic concepts of conservation of energy and momentum, causality, 
and Lorentz invariance that uniquely specifies the masses of the equivalent, bosons in the theory 
developed here. Despite the fact that the boson masses here are not constant and the values of 
the masses of the \\ and Z bosons can be quite different than those predicted by the Standard 
Model, tin* boson distribution functions are in perfect agreement with those in the parton theory 
for electroweak interactions. 


2.2. Bremsstrahlung Energy Loss 

Tin' calculation of the SP for bremsstrahlung processes involved the determination of several 
intermediate quantities. The most difficult one to derive was the frequency spectrum, (El/dAdu), 
of the field of bosons outside the charge; this quantity was referred to above as the boson distribu- 
tion function. < El/dAd A represents the total amount of energy carried by the field per unit boson 
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energy per area element, on the wavefront of the Lcrentz-c-emtrae-tecl E and D fields. In short, it is 
the Fourier transform of the Poynting ve'ctor S^ExD + M 2 $A mentioned above. 

A slight subtlety had to be addressed in deriving cl' 2 1 /dAdw. EM radiation is known to be 
polarized in a plane perpendicular to the direction of propagation, which means that when an EM 
wave (i.e., a swarm of photons) impinges on a charged particle, tin' particle will only move 1 in that 
transverse plain'. But radiation associated with massive' bosons can exist, in an additional state 
longitudinal polarization, which means that, the induced motion can take place m practically any 
direction relative to the direction of propagation of the wave. So there' are actually two different 
freciue'iiey spectrum functions --one corresponding to transverse'ly-polarize'd bosons (designated wit h 
a. T) and the' other describing lougitudinally-polarizeel bosons (designated with an L); it can be 
shown that only the T boson states contribute to the bremsstrahlung process. The' frequency 
spectrum for T boson states is 


d 2 1 
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= U 2 K'HO, 


where 1 


anel 


Io = 


2 (Jy + <Ia 
b 2 v' 2 


— const 


* = &WA/*4 — 


( 6 ) 

( 7 ) 

( 8 ) 


/v j is a moelifie'el Bessel function of the' see-onel kinel, of order 1, h is the impact parameter in the 
collision, and u = A E is the angular freeiuene y of the boson. 

Next, the radiation cross section (RCS), d\/dw, which conveys the probability for a boson 
with a specified energy to be' emitted from one charge and subseeiuently absorbed by another, 
was determined. It, is, in a crude sense, the? product of dPl/dAdw and another cross, section, the 
scattering cross sec tion (SCS), da/cin, all integrated ewer a solid angle of 4tt steradiaus and the 
(infinite) area of the wavefront, of Lorentz contracted E and D fields. The' SCS can be' thought of 
as the time average' of the probability for the incident charge to be' deflected into a certain specified 
direction during the collision with another charge'. For the EM case', this quantity is referred to 
as the Thomson scattering cross section (TSCS). A generalized version of it was derived here' from 
hist, principle's: 2 

(1(7 ((/incident) ((/target) 1- 
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((/incident ) 2 is Q 2 for the incident, particle 1 (e.g., the incident cosmic ray), (garget) Q“ i° r ^ 10 taI S ot 
particle (e.g., an atom in the space'craft wall), V is the velocity of the pulse of bremsstrahlung radi- 
ation t'mitte'd, and 6 is the scattering angle of the pulse of bremsstrahlung radiation omitted. The 


o 


formula n'ducvs to tho familiar Thomson formula when the charge's are set. to the electric charge' c. 
and M is set to zoro. 


After these two quantities determining d\/duj were identified, it was found that, unlike the 
traelitiemal RCS, the' gonuralizod RCS e'ould not be' a simple' product of the' two above 1 - me'utiejne'el 
terms. After ironing out this technicality, a final formula was founel for the gonoralizod RCS that 
immediate'lv yields a generalized version of the* SP: 


2 oc 

ffi 1 (//incident ) ((/target )*~ f n^_j-2(c\ -2Mb 

dw ~ 3 M'j J V V ‘ Ujr 


( 10 ) 


where b mm is the minimum impact parameter effective in the collision. A detaileel analysis shows 
that, le>r any of the three type's of eleetroweak interactions of intorost, the overall behavior of the 
RCS is dominated by the u — > 0 limit. Given that the energy A E of the boson that, mediates the 
interaction is related to u, ! via the equation A E = u>, it can be saiel that the nonrolativistie bosons 
(i.e*., those' with A E <C Ef, or x <G 1) contribute the most to the generation of bremsstrahlung. The* 
bosem mass sche'me' them iui[)lies that those' bosons with At Atf make' the* gieate’st e'outribution. 
As such, the choice of b lnlu for the* massive mediator eases shemlel be> exactly the same* as that made 
for the* photon case. It is easily worked e>ut that b , uin = l/2M f v is the* appropriate choice. The* 
ge'neralize'el SP cannot be cast into a simple compact formula and necessarily involves a two-fold 
integral: 


dE\ 


hrem 
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wlieir A is tho iiuinhor density ot atoms in the target material, and u? max is the maximum angular 
frequency (allowed by conservation of energy) of the mediating boson. 


In contrast to the above equation, the formula for the traditional (massless mediator) SP is a 
one-dimensional integral. Numerical methods were needed to solve the integrals. Numerical con- 
vergence of the integral expression for the RCS was assured by comparing the numerical result with 
the expected analytic result for the photon- mediated bremsstrahlung process. Tho same integra- 
tion step size was then used to calculate the RCS for W- and Z-boson- mediated interactions. As 
mentioned above, the SP was then determined by performing one more integral, using the same 
integration step size as before. For ultrarelativist ie cosmic rays (in the complete screening limit), 
the RCS turned out to simply be a constant (independent of the energy of the cosmic ray) for the 
photon case, and very nearly a constant for the W- and Z-boson cases. As such, tho SP is the 
product of the RCS, the density of the bulk matter, and the energy of the cosmic ray. 

Written in another more intuitive way, the energy E( X) of the cosmic ray can be shown to be 
related to the distance traversed X as: E(X) = Here, Eq is the energy of the cosmic 

ray as it enters the material and X {) is a parameter known as the radiation length. Tho radiation 


G 



length is the distance it takes for the particle’s energy to fall to 1/e of its initial value. According 
to an oversimplified analysis (found in the literature) of EM brenisstrahlung emitted by an electron 
traversing three different possible materials (lead, aluminum, and air), X Q is about 0.4 cm for lead, 
7 cm for aluminum, and 270 m for air. These results were verified, and corresponding valla's of X„ 
for both W- and Z-boson-mediated processes were determined. For lead, A’ () was found to be 45 
m and 7.9 cm for the W- and Z-boson cases, respectively; for aluminum, A',, was found to be 20 
m and 3.3 m for the W- and Z-boson cases, respectively; and for air, A« was found to be 20 km 
and 14 km for the W- and Z-boson cases, respectively. As can be gleaned from the above formula, 
smaller values of A„ correspond to greater amounts of energy lost in the form of brenisstrahlung. 
So clearly, the greatest amount of brenisstrahlung energy is carried away by photons and the least, 
amount is carried away by W-bosons. In fact, the proportions of brenisstrahlung .energies carried 
away by these three particles can be easily determined. If an electron is the incident, cosmic ray 
(traveling at any speed), the photon : Z- boson : W-bosou ratios of bremsstrahlung energy loss are 
roughly 1 : 10“ 2 : 10 r ’ for lead, 1 : 10“ 2 : 10~ 4 for aluminum, and 1 : 10' 2 : HE 2 for air. 

The cases considered in this study assume the cosmic ray to be an electron because electrons 
are the most efficient, of all cosmic rays in emitting brenisstrahlung. As a rule of thumb, greater 
amounts of brenisstrahlung will bo produced by cosmic rays with higher charge-to-mass ratios 
traveling through materials with greater charges. The most efficient scenario turns out to be an 
electron traveling through a dense material such as lead. In figure 1, a comparison is made between 
the fraction of cosmic ray energy lost and brenisstrahlung associated with the three elect loweak 
bosons for a 500 GeV electron traveling through a 0.4 cm slab of lead. As alluded to above, the 
energy lost, to W- and Z-boson bremsstrahlung is a mere 0.01 and 4 percent, respectively, of that lost 
to EM bremsstrahlung. So, it can be reasonably stated that the overwhelming amount of energy 
lost in the form of brenisstrahlung is carried away by photons. 

2.3. Collisional Energy Loss 

The calculation of the SP for collisional energy loss processes was easier than that for 
bremsstrahlung processes. Standard Methods (in particular, Fourier analysis) were employed to 
find the amount of energy E transferred to an electron orbiting an atom in a grazing collision wit h 
a fast- moving (incident) charge, 


E = 


incident ) ((/electron) 


1 {<] 

8tT 2 M, lech,, n l’ 2 ^ 2 




A pu 
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where A p ±min = v/yv, A p ± = yj A I 2 + (Apx mi „) 2 , and £ = bAp±. It was a simple matter to 
convert this generalized E into a generalized SP. The SP found for these collisional processes was 
much simpler in form than that for bremsstrahlung processes; unlike the bremsstrahlung case, the 
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SP hero reduced to a simple compact formula, albeit, quite complicated: 
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where 0 — ^niinA/>j_. When Q is set to the electric charge and AI is set to zero, the formula reduces 
to the familiar traditional version. 


One of the main goals of this paper is to compare the amount of energy lost via this mechanism 
with that lost via bremsstrahlung for the three different types of electroweak interaction. It is 
known that the former energy loss mechanism dominates t.lie latter for slower moving electrons. 
Again, an electron is taken to be the incident cosmic ray particle because it is the most efficient 
of all particles at radiating energy. Specifically, it is known that collisional energy loss dominates 
bremsstrahlung energy loss if the relativistic factor 7 (energy-to-mass ratio) of the cosmic ray is 
less than about 20 if the bulk matter through which the electron travels is lead, and 200 if the 
material is air. These thresholds (for EM processes) were verified in the present analysis, though 
the values found were about 17 and 170, respectively, and corresponding ones were found for W- 
and Z-boson processes. I11 lead, the threshold value of 7 was found to be about 8.4 x 10 1 and 17 
for \\ - and Z-boson processes, respectively. In air, the threshold value of 7 was found to be about 
5 x h) * and 1 x 1(P for W- and Z-boson processes, respectively. The results, assuming the cosmic 
ray to be an electron and the target to be made of lead, are shown in figure 2, where the ratio of 
the SP for the bremsstrahlung process over that for collisional energy loss process is plotted against 
7. Tin* threshold above wherein bremsstrahlung dominates collisional energy loss corresponds to 1 
along the y-axis. This threshold is clearly surpassed for all three types of electroweak interactions 
above collision energies corresponding to 7 = Iff’. The conclusion to be drawn from these results, 
then, is the fact that for ultrarelativistic cosmic rays, where 7 is typically greater than about 1(P. 
the dominant energy loss mechanism is most assuredly bremsstrahlung, instead of collisions, for all 
three types of electroweak interactions. 

3. Conclusions 

In summary, the stopping power for a cosmic ray traveling through bulk matter was derived tak- 
ing into account both collisional and Bremsstrahlung energy-loss mechanisms. Numerical routines 
were used to evaluate the stopping power for different scenarios. In all cases, the cosmic ray was 
assumed to be an electron, and three different possible materials were considered as the bulk mat- 
tei. lead, aluminum, and air. One major result of the project was that bremsstrahlung dominates 
the collisional energy- loss mechanism for the ultrarelativistic cosmic rays of interest for all three 
different types of electroweak interactions. The other major result showed that the bremsstrahlung 
associated with W- and Z-bosons only enters as a correction at the 1-percent level to electromag- 
netic bremsstrahlung. This latter result is the more important, one for the purposes of figuring out 
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the influence that the most, energetic forms of space radiation have on humans. It shows that, weak 
force effects can be safely ignored compared with similar electromagnetic effects. 
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